Spectral signatures of the Fulde-Ferrell-Larkin-Ovchinnikov order parameter in 

one-dimensional optical lattices 
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We address an imbalanced two-component atomic Fermi gas restricted by a one-dimensional (ID) 
optical lattice and an external harmonic potential, within the mean-field Bogoliubov-de Gennes 
formalism. We show that characteristic features of the Fulde-Ferrell-Larkin-Ovchirmikov state are 
visible in the RF-spectra and in the momentum resolved photoemission spectra of the gas. Spe- 
cially, Andreev states or mid-gap states can be clearly resolved, which gives a direct experimentally 
observable signature of the oscillating order parameter. 
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Experimental realization of spin-density imbalanced 
Fermi gases [H, d, Q has opened exciting new possibil- 
ities to study pairing in systems where matching of spin- 
resolved Fermi surfaces, which is the base of Bardeen- 
Cooper-Shreffer (BCS) theory, is not valid. Either phase 
separation or extension of the BCS pairing to some other, 
exotic mechanism is inevitable. The question is of inter- 
est in context of various solid state materials 0, [H, [f| 0] as 
well as in hight-energy and astrophysics Q. One of the 
main candidates for the non-BCS pairing is the so-called 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state @, [l(J. 
The zero temperature ground state properties of the 
FFLO state in the context of ID Fermi gases have been 
studied extensively within mean-field, exact, and DMRG 
approaches The rapid development of exp eriments 
on Fermi gases in optical lattices 12, [3, H EH sug- 
gests that such systems will be available soon. How- 
ever, it is not obvious how to observe e.g. the spatially 
varying order parameter and other characteristics of the 
FFLO state, although noise correlations have been pro- 
posed to provide information about the pairing correla- 
tions [16]. In this letter we show how the characteris- 
tics of the FFLO state are prominently reflected in RF- 
spectroscopy 17] and in the recently introduced [181 ] pho- 
toemission spectroscopy of Fermi gases. Unlike the BCS, 
the FFLO state allows population of single-particle ex- 
citations even at zero temperature, corresponding to un- 
paired particles. For a spatially non-uniform order pa- 
rameter, some of these excitations can be understood 
as Andreev bound states residing close to the nodes of 
the order parameter. We show that such excitations pro- 
duce distinct features in the spectra, at negative energies, 
and thus provide a signature of oscillations of the order 
parameter that is easily distinguishable from the usual 
pairing signatures at positive energies. Furthermore, by 
calculating spectra also at finite temperatures, we show 
that such features are uniquely related to oscillations of 
the order parameter. 

We consider a two-component attractive Fermi gas 
confined by an external potential in a ID lattice with 



L sites. At low filling, this corresponds to a one- 
dimensional gas without the lattice. Our qualitative re- 
sults concerning the spectral signatures should be valid 
in this case as well. We apply a mean-field Bogoliubov- 
de Gennes (BdG) approach. In ID, long range order is 
absent and the mean-field approximation is not as well 
valid as in 3D, especially the value of the mean- field order 
parameter may deviate from the exact value [l9j . How- 
ever, mean-field approaches are often applied since they 
provide qualitative information on the system, as well as 
make it possible to study the effect of finite temperature 
and to take into account non-trivial confining geometries 
(which is not possible in the case of exact analytical solu- 
tions) . Here, the mean-field description further allows to 
calculate the RF-spectrum as well as to understand and 
interpret it in a transparent way. 

At the mean-field level, the system is 
by the single-band Hubbard-Hamiltonian: 

(4 + h - c ) + £i(A4r4i ■ 

J2ia ((^ cxt ~ AV) - Ufiis^ h ia . Here the spin labels a =\ 
, J. refer to two hyperfine states of the atoms, the pairing 
gap is defined as A, = — [/(cjj, Cjj), where U{> 0) de- 
notes the on-site attractive interaction, t is the hopping 
strength and I^ oxt = Vo(i — L/2) 2 is an external harmonic 
potential. The a refers to opposite spin component of a. 
We neglect the Hartree interaction term Ufii S because 
doing so does not significantly affect the results but sim- 
plifies the numerical calculation. We use the Bogoliubov 
transformation, c irJ = J2 a ( u aicr Tew ~ w* aicr tf aS ) to diag- 
onalizc the Hamiltonian TL, leading to the BdG equation 
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gle particle Hamiltonian and Aij = Aj<5;j de- 
notes the local pairing gap. The above equation is 
solved together with the self-consistency conditions 
for the gap and the spin-resolved densities: nf — 
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r \ 2 n F (E aa ) + \v aia \ 2 n F (-E aS ) 



u a i]V ai in F (E a ^) - u a ii v a ^n F (-E al ) 

Here the densities are normalized to iVf and N±, and n F 
is the Fermi function at the temperature T. 

In RF spectroscopy, the internal state of one of the 
components (f or j, energy w vi ) is coupled to a third 
internal state (denoted final state, energy ujf) by an RF 
pulse. The number of particles transferred from the ini- 
tial to the (initially empty) final state can be observed 
and gives the spectrum as a function of the detuning, 
J UI (S), where S = lurf — (u>f — w T/i ). Recently, first 
experiments where the final state momentum was re- 
solved were performed [l8j |. in analogy to photoemission 
spectroscopy of solid state systems. The corresponding 
spectrum is then J rn (S, k) where k is the final state mo- 
mentum. We have derived J r/l (S,k) and J T / , (5) in the 
present case, following the approach [20L |21|. This cor- 
responds to the quasiparticle picture [22j where the RF 
field can be understood to create quasiparticle excitations 
by pair breaking rather than to induce coherent rotation 
of the internal states. It is valid when the final state in- 
teractions are negligible [23|, 0]. The spectra become: 

J-t/i( S ) =ELi J t/i(^) and 

L N L 2 

a— 1 i=X 

x S(E ain + e K - 6 - jti m ) 
L 2 

(Earn) 



x S(E aVl - e K + 5 + fi rn ) 



(1) 



Here the symbol f / J. denotes the spin-resolved case, 
while jj. and "non" refer to the balanced, non-interacting 
case. The particles in the final state are assumed to be 
trapped in the same potential, therefore the final state 
wavefunctions in the overlap integrals, v™*i,, are calcu- 
lated by solving the BdG equations in the non-interacting 
case, for high filling in order to have all relevant wave- 
functions v™^ non-zero. The final state dispersion is also 
obtained from such calculation: ejc = /x n — E Klv In the 
numerical calculations, we replace the <5-function by a 
Lorentzian distribution C{x) = ^ x 'i S j rF 'i an d use the line 
width r = 0.05. All energies are in units of t, and we 
use a realistic trapping potential Vo = 5 x 10~ 4 and the 
attractive interaction strength U = —3. 

Note that in Eq.(p} K is the quantum number of the 
eigenstate of the combined lattice plus harmonic po- 
tential. In the limit of vanishing harmonic potential, 
it approaches the lattice momentum k. In an experi- 
ment, it is the lattice momentum that can be resolved 
from absorption images after the time of flight expan- 
sion. Therefore we calculate the observable momentum- 
resolved spectrum by taking into account the overlap 
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FIG. 1: (color online) Zero-temperature densities and order 
parameters A as a function of site position i for polarization 
P = 0.23 (iV T = 40, Ni = 25) (top panel) and P = 0.7 (40, 7) 
(bottom panel) in L = 150 lattice sites, and in the presence 
of a harmonic potential. Top panel shows A also at finite 
temperatures Tl — 0.075 and T2 = 0.1. The wavefunction 
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is also shown in bottom panel. 



between the eigenfunctions for the quantum number K 
and those corresponding to the lattice momentum k: 

J VI {8,K). Effec- 



tively, this describes the averaging by the external har- 
monic potential of the momentum resolved spectra. 

Typical density and order parameter profiles are shown 
in Fig. [Tj for the polarizations P = 0.23 (iV T = 40, AT X = 
25) and P — 0.70 (40, 7) where the polarization is defined 
as P = (7V T - iV x )/(iV T + Ni). For P = 0.23, the order 
parameter is constant and the density imbalance zero in 
the center of trap. Surrounding that, there is an area of 
oscillating order parameter and finite density imbalance. 
For P = 0.7, this area spans through the central part 
as well. At the edges, in both cases, the order param- 
eter vanishes while the majority particle density is still 
non-zero. This kind of BCS-FFLO-Normal and FFLO- 
Normal shell structures are in accordance with ll|. For 
the parameters used here, the BCS part vanishes and 
the FFLO area reaches the trap center for polarizations 
P > 0.3. 

Figure [2] shows the RF spectra for different polariza- 
tions. Both the minority and majority spectra have a 
peak located at the same, positive, detuning. This cor- 
responds to breaking of pairs, in other words, excita- 
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FIG. 2: (color online) Majority and minority (inset) RF spec- 
tra for different polarizations. The curve marked with open 
circles is explained in the text. 



tions are created in the gas. The majority spectrum 
displays also a prominent peak at zero detuning, cor- 
responding to unpaired majority component particles. 
Note that there is considerable spectral weight at neg- 
ative detunings. We also plot separately the second term 
(proportional to n-p(E al/l )) in Eq.|T]) for the polariza- 
tion P = 0.23 (open circle): this gives almost all of the 
spectral signatures at negative detunings, and part of 
the zero detuning peak for which also the first term con- 
tributes. In both terms, the energy conservation relation 
imposes 5 — \Ei a \ — — ck)- For a = K this is al- 
ways positive or zero, for the normal and BCS states. 
The BdG wavefunctions v and u in the interacting case 
are, however, not the same as in the non-interacting case 
and therefore the overlap integrals in Eq. ((]} can be non- 
zero also for a ^ K. Thus, in order to obtain spec- 
tral weight at negative detunings, the interacting case 
wavefunctions have to be decomposed of a large set of 
non-interacting case wavefunctions. This may happen 
when the wavefunctions have to accommodate to a spa- 
tially highly non-uniform order parameter. We have an- 
alyzed the energy and eigenvalue structure of the sys- 
tem and found that those states that contribute to the 
negative detuning peaks correspond to majority particle 
wavefunctions that have maxima at the nodes of the or- 
der parameter. As examples, we discuss here the states 
a = 142 in the P = 0.23 case, and a = 123 for P = 0.7. 
In both cases, these states correspond to negative but 
close to zero excitation energies, i.e. they are close to the 
Fermi surfaces. By calculating the overlap integrals, we 
found that the state a = 142 (P = 0.23) spans all the way 
from K = 100 to K = 140. From the energy conservation 
relation one can then see that e.g. the K = 100 contribu- 
tion leads to spectral weight at 5 ~ —0.6, and the 140 at 
5 ~ —0.3. Similarly, for P — 0.7, the state a = 123 is de- 
composed of K from 70 to 150, leading to contributions 
as far as 5 ~ — 1. The state a = 123 is shown in Figure 
1. It resides inside the order parameter, with maxima 
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FIG. 3: (color online) Final state momentum resolved spectra 
for majority and minority components for different k (top 
panel), for P — 0.7. Majority component spectra at different 
temperatures for P — 0.23 (bottom panel). 



at the nodes. Typically, states contributing to negative 
and zero detuning were found to be located either totally 
inside the order parameter, or partly outside and partly 
inside, and always following the node structure. 

Some of the states discussed above are related to An- 
dreev bound states formed at FFLO superconductor - 
ferromagnet interfaces [25j. In semiclassical thinking, 
the spatially varying order parameter forms a potential 
affecting the unpaired majority particle excitation spec- 
trum. In solid state systems, the Andreev bound states 
often lead to zero bias anomalies in tunneling experi- 
ments. RF-spectroscopy, in the quasiparticle picture as 
considered here, is analogous to tunneling and zero bias 
corresponds to zero detuning. However, part of the spec- 
tral signatures of the Andreev type states are at neg- 
ative detunings due to the momentum conservation in 
the spectroscopy. Without momentum conservation, the 
overlap integrals loose their significance and transitions 
to all final states are allowed, then the densities of state 
of the final and initial state will determine the spectral 
weight. Since density of states typically grows for higher 
lying states, this would favour zero detuning transitions 
since the negative detuning transitions are to lower lying 
final states. The lower lying final states are often Pauli 
blocked in tunneling experiments. 

In Fig. [31 we show also the fc-resolved spectra that can 
be observed by photoemission spectroscopy. The final 
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states around the Fermi level of the majority component 
correspond mainly to the non-paired particles, whereas 
the low momenta are involved in both pairing and non- 
paired excitations. Note that the momentum-resolved 
RF-spectroscopy may help to distinguish negative detun- 
ing contributions since peaks instead of a steady pedestal 
appear. 

It is now important to discuss the dependence of these 
negative detuning spectral features on the system size. 
We have varied the lattice size and particle number from 
L = 100 -> 150, N = 40 — > 80 and the results are similar; 
from this and from the understanding of the results given 
above it is obvious that they are not any numerical arti- 
fact related to a finite size system. They are, however, a 
mesoscopic effect in the sense that the order parameter 
oscillations have to be in a scale that, in the semiclassical 
thinking, provides sufficiently tight confinement to obtain 
Andreev states whose spread in eigenstates of the poten- 
tial is large. The spread has to be large enough to allow 
transitions where a and K states are sufficiently different 
in energy compared to other broadenings in the exper- 
iment. Note that in 2^, 27 1 it was shown that a peak 
at zero detuning in the minority spectrum is a signature 
of order parameter oscillations; this applies to large sys- 
tems and for high polarizations, and is not visible in the 
present case. Ultracold atoms in optical lattices provide 
a system where ID potentials accommodating, e.g., a few 
hundreds or thousands of particles can be realized. Such 
numbers are not essentially different from those used in 
our simulations, therefore we believe that spectral weight 
on the negative detunings should be visible in the exper- 
iment. 

Finally, although it is clear from the above discussion 
that the negative detuning features require the existence 
of a spatially non-uniform order parameter, one might 
still ask whether oscillations are required or is simply a 
smoothly varying order parameter profile sufficient? To 
clarify this, we calculated spectra also at finite temper- 
atures. The disappearance of the order parameter os- 
cillations with finite temperature is shown in FigQ] and 
the finite T spectra in FigGH for P = 0.23. The pedestal- 
type spectrum at negative detunings disappears when the 
oscillations vanish. Also the zero detuning peak eventu- 
ally goes down since finite temperature enables pairing 
even in the presence of population imbalance, i.e. the gas 
enters the finite temperature BCS phase. We can thus 
conclude that the prominent, pedestal- like spectrum at 
negative detunings, with a peak at zero, is a unique sig- 
nature of an oscillating order parameter. This would be 
one of the most direct signatures of the FFLO state ob- 
servable in any system so far. 

In summary, we considered the density imbalanced, at- 
tractively interacting Fermi gas in a ID optical lattice, 
using the mean-field BdG formalism. We propose that in 
such systems, RF-spectroscopy and photoemission spec- 
troscopy can provide information not only about the pair- 



ing gap and the amount of paired/unpaired minority and 
majority particles, but also about the spatial structure of 
the order parameter. This is basically due to the momen- 
tum conservation in the process which gives an important 
role to wavefunction overlaps in transition probabilities 
and thus allows to resolve also spatial features of exci- 
tations. By analyzing spectra and eigenstates both at 
zero and finite temperatures, we have shown that sig- 
nificant spectral weight on negative detunings is a di- 
rect signature of the Andreev bound states or mid-gap 
states formed due to a spatially strongly modulated order 
parameter, and thereby a direct signature of the FFLO 
state. 
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